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CD ■ Abstract 



We discuss the effects of one additional sterile neutrino at the Neutrino Factory. Compared 



>• ' to earlier analyses, which have been motivated by LSND results, we do not impose any 
^2 • constraint on the additional mass squared splitting. This means that the additional mass 
"Tj- ■ eigenstate could, with small mixings, be located among the known ones, as it is suggested by 
^ ! the recent analysis of cosmological data. We use a self-consistent framework at the Neutrino 
^ I Factory without any constraints on the new parameters. We demonstrate for a combined 

(^ I short and long baseline setup that near detectors can provide the expected sensitivity at 
the LSND-motivated A?724]^-range, while some sensitivity can also be obtained in the region 
of the atmospheric mass splitting from the long baselines. We point out that limits on 
K^ ! such very light sterile neutrinos may also be obtained from a re-analysis of atmospheric 

H [ and solar neutrino oscillation data, as well as from supernova neutrino observations. In 
the second part of the analysis, we compare our sensitivity with the existing literature 
using additional assumptions, such as |A?n4^| ^ |Am|;^| leading to averaging of the fast 
oscillations in the far detectors. We demonstrate that while the Neutrino Factory has 
excellent sensitivity compared to existing studies using similar assumptions, one has to be 
very careful interpreting these results for a combined short and long baseline setup where 
oscillations could occur in the near detectors. We also test the impact of additional Ur 
detectors at the short and long baselines, and we do not find a substantial improvement of 
the sensitivities. 
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1 Introduction 

Neutrino oscillation experiments recently have given us compelling evidence that active neu- 
trinos are massive particles [1], pointing towards physics beyond the Standard Model. In 
the standard three-generation scenario, there are two characteristic mass squared splittings 
{Am'^i , Am^i) and three mixing angles {612 ,^13 ,^23) as well as a CP violation phase ^cp 
affecting neutrino oscillations. Disappearance of muon neutrinos, which is mainly driven 
by lArrig^l and ^23? has been observed in the atmospheric neutrino oscillation experiments, 
such as Super-Kamiokande [2], and in the MINOS long baseline experiment [3]. Disappear- 
ance of electron neutrinos has been observed from solar neutrino oscillation experiments 
very sensitive to 612 [4], whereas A??7,^2 has been strongly constrained by the KamLAND 
long baseline reactor neutrino experiment [5]. The CHOOZ short-baseline reactor neutrino 
oscillation experiment [6] has provided a limit sin^ 2^13 < 0.1. There are still unknown 
questions in the standard scenario: Am|]^ > (normal ordering) or Am|^ < (inverted 
ordering); the value of ^13, as there has been a recent hint for 6*13 > [7], and whether there 
is CP violation (CPV) in the lepton sector. Apart from these measurements, there has 
been the exceptional LSND measurement with an incompatible anomaly [8]. The simplest 
interpretation has been an additional sterile neutrino added to the standard picture with 
|Am4^| ^ |Am|^|. A global fit to all experimental data, however, is not in favor to this 
hypothesis [9], which means that more exotic scenarios would be required to describe this 
anomaly, such as a decaying sterile neutrino [10]. The recent results from MiniBooNE, how- 
ever, are consistent with sterile neutrino oscillations in the antineutrino sector [11]. Note 
that the LSND interpretation requires significant mixings with the active neutrinos, whereas 
small ad-mixtures, even if lAm^]^! ^ |A?Ti3^|, are not excluded. On the other hand, sterile 
neutrinos with IAtti^^I ~ lATTig-,^! or lAm^-,^! ~ Am^i, as they are motivated by a recent 
cosmological data analysis [12], have been hardly studied in the literature. Of course, such 
sterile neutrinos have to have small mixings with the active ones in order not to spoil the 
leading three-flavor fits. Therefore, we discuss sterile neutrinos beyond LSND, i.e., with- 
out any constraints to Aml^ including the full range. We consider the simplest case of 
only one additional sterile neutrino. Note that we do not impose cosmological constraints 
on the sterile neutrino masses, but we assume that the masses are to be constrained in a 
self-consistent way with terrestrial experiments. 

Since we cannot rely on constraints from the atmospheric and solar experiments, which 
need to be re-analyzed in a global fit for the presence of very light sterile neutrinos, we 
consider the Neutrino Factory (NF). This instrument provides excellent sensitivities to the 
standard oscillation parameters, and it can be used for a self-consistent simulation. In 
the proposed experiment, neutrinos are produced from muon decays in straight sections 
of a storage ring. The design of the Neutrino Factory has been extensively discussed in 
international studies, such as in Refs. [13-15]. The International Neutrino Factory and 
Superbeam Scoping Study [15-17] has laid the foundations for the currently ongoing Design 
Study for the Neutrino Factory (IDS-NF) [18]. This initiative from about 2007 to 2013 is 
to present a design report, schedule, cost estimate, and risk assessment for the Neutrino 
Factory. IDS-NF defines a baseline setup of a high energy neutrino factory with E^^ = 
25 GeV and two baselines Li ~ 4 000 km and L2 — 7500 km (the "magic" baseline [19]) 



operated by two racetrack-shaped storage rings, where the muon energy is 25 GeV (for 
optimization questions, see, e.g., Refs. [20-27]). A key component is the magnetized iron 
detector (MIND) as far detector, where the magnetization is necessary to distinguish the 
"right-sign" {e.g., from u^ — )■ i/^) from the "wrong-sign" {e.g., from z/g ~^ '^fi) rnuons. 
The identification of the muon charge of the wrong-sign muons allows for CP violation 
measurements in the muon neutrino appearance channels [21,22]. Physics with near detector 
configurations for cross section and flux measurements have been discussed in Refs. [16, 
28]. Note that recently a smaller scale low energy version of a neutrino factory is also 
been proposed with E'^ ~ 5 GeV in Refs. [29-34]. The Neutrino Factory is claimed to be 
a precision instrument not only because it can answer the unknown questions addressed 
above, but also because it can tell us the story beyond three flavor neutrino oscillation 
physics. Examples are unitarity violation of the mixing matrix coming from heavy fermion 
singlets [35-40] and non-standard interactions during neutrino productions, propagation in 
matter, or detection [27,40,41,43-47]. In this study, we focus on the constraint of light 
sterile neutrinos at the Neutrino Factory, for earlier studies, see Refs. [38,48-55]. These 
are either based on the short-baseline (large Am^^ dominance) or long baseline (large Amlj^ 
effects average out) limit. We combine these approaches in one analysis, including possible 
near detectors and additional channels. Note that in certain short-baseline channels CPT 
invariance tests are possible [54], as they may be motivated by the recent MiniBooNE 
results [11]. 

The structure of the paper is organized as follows. In Sec. 2 we illustrate and discuss the 
parametrization of the neutrino mixing matrix and show approximate expressions for four 
generation neutrino oscillation probabilities. Matter effects are taken into account in the 
long baseline setup. We also discuss the experiment simulation. In Sec. 3, we show the most 
general limits on the new parameters without any additional assumptions. A comparison 
with earlier studies using specific assumptions and a discussion of related issues can be 
found in Sec. 4. Finally, we summarize in Sec. 5. 

2 Analytical discussion and simulation techniques 

In this section, we first discuss the parametrization of the mixing matrix. Then we de- 
rive analytical formulas useful for a first qualitative discussion. Finally, we describe our 
simulation techniques and assumptions. 

2.1 Mass schemes and parametrization of the mixing matrix 

We study the constraints to the general four neutrino scheme without any assumptions at 
the Neutrino Factory. While the global fits imply the mixings between active and sterile 
neutrinos must be small, the four neutrino scheme has to recover the standard picture once 
we switch off the small mixings between active and sterile neutrinos. The four neutrino 
schemes can be categorized into two different classes: the 2+2 scheme, in which the solar 
and atmospheric mass squared splittings are separated by a new splitting, and the 3+1 
scheme, in which the new mass eigenstate is added somewhere to the existing mass pattern. 
The 2+2 scheme is, at least for an LSND-like new mass squared splitting, strongly disfavored 
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Figure 1: The mass ordering of four neutrino eigenstates (not to scale). The arrow ihustrates the 
new characteristic mass squared difference Aml^. The four different scenarios correspond to Ayrigj^ > 0, 
Amli > (A), Am§i > 0, Ato|i < (B), Am§i < 0, Am|i > (C), and Am§i > 0, Am^i > (D). 



by global fits [56,57]. The 3+1 scheme, on the other hand, naturally recovers the standard 
picture in the case of small mixings. Therefore, naturally, we consider the 3+1 scheme 
only. We show the possible mass ordering of four neutrino eigenstates in Fig. 1. The arrow 
illustrates the new characteristic mass squared difference Amii. The four different scenarios 
correspond to Amj^ > 0, Am^^ > (A), Am^^ > 0, Amj^ < (B), Am^^ < 0, Am^^ > 
(C), and Am|]^ > 0, Amli > (D). Unless noted explicitly, we show the results for scenario 
(A). 

As far as the parametrization of the four flavor mixing matrix Uai is concerned, the straight- 
forward choice is a parametrization-independent approach. For example, some of the oscilla- 
tion probabilities in the short-baseline (one mass squared dominance) limit IA41I ~ 0{1) ^ 
IA31I with Aij = Am|L/(4E) read [55] 
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Information on the combination 4|[/^4p|?7^4p, for instance, can be easily obtained from the 
NOMAD [58] and CHORUS [59] experiments. The mixing matrix elements U^4 and f/^4 
can then be interpreted in any arbitrary parametrization of U. Suppose now we wanted 
to constrain each mixing matrix element |t/ajp individually. It is obvious from the above 
equations that |t/^4p could be constrained by V^^. The same applies to |?7e4p by Vee- On 
the other hand, let's take a look at |t/T-4p- It would be directly accessible by Vrr, which is, 
of course, basically inaccessible by any experiment. The alternatives are Ver and V^r, where 
|[/^4p enters as a product with |f/e4p and |f/;i4p, respectively. If one simulates |t/r4p = 
to obtain the expected bound for |t/r4p! the marginalization over these two other elements, 
ending up at zero, will inevitably destroy the sensitivity. Therefore, unless |f/e4p > or 
1^^*41^ > in Nature, one cannot obtain a bound on the individual |t/T4p from these short- 
baseline probabilities. In summary, one cannot constrain all of the individual elements 
at the short baseline. We will recover this parametrization-independent statement in our 
parametrization below. 



Another comphcation of the parametrization-independent approach is that it does not allow 
for the simultaneous treatment of short and long baselines, and for arbitrary values of Amli. 
For example, the mixing matrix element combinations in the above equations obviously 
come from the CP conserving part of the probabilities and rely on the fact that one mass 
squared splitting is much larger than the other two. If arbitrary Amli are allowed, other 
combinations of mixing matrix elements enter, and, in principle, the individual elements 
are needed separately even in the absence of CP violation. Furthermore, charged and 
neutral current matter effects at the long baselines complicate this picture even further, see 
Ref. [55] . Choosing a particular parametrization for the unitary mixing matrix is therefore 
the commonly used approach, even in the three flavor picture. Compared to using the 
mixing matrix elements individually, it has the advantage that unitarity is automatically 
conserved. Of course, the results are then parametrization-dependent, but the qualitative 
features are the same as in the parametrization-independent approach, as we will show 
below. 

In principle, there are many different parametrization of the neutrino mixing matrix, since 
the order of the sub-rotations is arbitrary. For the sake of simplicity, we impose the following 
requirements: 

• The standard PMNS matrix has to be recovered in the case of small new mixing 
angles; that fixes the order of the corresponding sub-sector rotations 

• The phases are attached to the 12-, 13- and 23-rotations. Therefore, if one of the 
standard mixing angles can be rotated away in a particular measurement, the corre- 
sponding phase also automatically becomes unphysical. 

• The order of the 34-24- 14-rotations is arbitrary. We choose the 34-angle as the left- 
most one, which makes it hardest to observe (it affects the i/7--i/s-mixing). Changing 
the order here does not change the fact that one of the rotations is difficult to extract. 

Our parametrization explicitly reads 

U = i?34(^34, 0) i?24(^24, 0) i?i4(^14, 0) i?23(^23, ^s) i?13(^13, S^) RuiOu, ^i) . (5) 

In Eq. (5), Rij{6ij, 6i) are the complex rotation matrices in the zj'-plane, defined as: 
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This means that ^2 becomes Sqp in the three fiavor limit. Note that the order between 
the 14-and 23-rotations is arbitrary, since these matrices commute. We do not involve any 
Majorana phases for their absence in the oscillation probabilities. The parametrization in 
Eq. (5) is the same as the one used in Refs. [53,60], and the same as the one in Ref. [61], 
since they fix their 62, corresponding to our 61, to zero. 



2.2 Analytical considerations 

The general evolution equation of flavor eigenstates in matter can be expressed as 

^j^W) = n^^Wp) = ^{UD'U^ + A)\up) , (7) 

where D = Diag{mi, m2, ma, 1714} is the diagonal mass matrix and A = Diagjoe, 0, 0, — a„} 
is the matter potential with Ue = 2ylEGFne and a„, = \/2EG pfin (see, e.g.^ Ref. [62] for 
details). Here Oe and a„ are the electron and neutron number density in matter, respectively, 
with Ue — Un in Earth matter. Note that the neutral current matter effect may lead to 
additional matter-driven effects, which are not easy to capture in an analytical treatment. 
The perturbative solution of this equation depends on the A4i-regime. We will discuss 
different regimes separately. 

In the one mass squared dominance (short-baseline) limit IA41I ~ 0{1) ^ IA31I, the matter 
effects can be safely ignored, and the probabilities read: 



(9) 
(10) 

(11) 
(12) 

V^.s = 4C14C24C34S24 sin^ A41 (13) 

^MM = 1 - cf^sl^ [3 + 24 COS (2^24) - COS (2^14)] sin^ A41 (14) 

Note that the neutral current rate is proportional to 1 — Ves for electron neutrinos at the 
source, and 1 — P^s for muon neutrinos at the source. These probabilities correspond to 
the near detector limit at the Neutrino Factory, where the baseline depends on the sensitive 
|Am4^|. For example, at (i ^ 18 km, one has optimal sensitivity for |A?ti|j^| ~ leV^, whereas 
at (i ~ 1.5 km, one has optimal sensitivity for |A?7i4^| ~ 10 eV^ with the 25 GeV Neutrino 
Factory (c/.. Fig. 2 in Ref. [54]). Here d is the distance to the end of the decay straight, 
which is related to an effective baseline L^s = \/d{d + s) with s the length of the decay 
straight if the decays are averaged over the straight [28]. 

Let us now compare these equations with the parametrization independent approach in 
Sec. 2.1, Eqs. (1) to (4), assuming that the new mixing angles are small. Obviously, we 
recover the same qualitative features: ^24 can be measured by V^^ in the short-baseline 
limit, and ^14 by Vee- On the other hand, ^34 only shows up in combination with the 
other small mixing angles. This means that by the same arguments as in Sec. 2.1 in the 
parametrization-independent approach for |t/T4p; no separate limit on ^34 can be obtained 
at the short baselines. In the long-baseline discussion, this parameter is therefore the most 
interesting one. Note that apart from the above one mass squared dominance (|A?Ti4]^|) 
discussion, also the two mass squared dominance (|Am|]^| and |A??7,|^|) at somewhat longer 
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baselines (L > 10 km at E^ = 25 GeV) is interesting, because CP violation effects may be 
observable there [51]. However, we focus on near detectors closer to the source in this study. 

For long baselines, some of the relevant features of the probability transitions can be well 
understood using simple perturbative expansions in appropriate small dimensionless param- 
eters with A31 = C(l) ^ A41, i.e., if the effects of the large |Am4;^| are averaged out. In 
particular, up to the second order in S13, S14, S24, S34, ,§23 = sin 6*23 — 75, and considering A21 



as small as s?-, we obtain ^: 
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^Ms =2-^24(3 + cos(2A3i)) + S34sin^(A3i) - S24S34 sin(53 sin(2A3i) , (20) 

where we used the short-hand notation Ae„ = a^^nL / {AE) . From these formulas we learn 
that at the long baselines, ^14 is difficult to measure, because it shows up to leading order 
only in Vee as constant term, which we do not consider at the long baseline. On the other 
hand, ^24 is best accessible by V^^ with the first term proportional to cos^(Ai3). From the 
above equations, the leading sensitivity to 6'34 can be expected from V^t at a long baseline, 
as claimed by Ref. [53]. However, the statistics in this channel is very low, which means 
that higher order terms in the other channels could still provide the leading sensitivity. 
Depending on the performance indicators, one of the problems is also the simultaneous 
presence of ^24 and S23 with the same spectral dependence in the first line of "P^^- Thus, 
if ^24 is marginalized over, the sensitivity to 6'34 disappears. We will test these effects 
quantitatively later. 

Another interesting case, which can be treated analytically, is IA41I ^ A31, i.e., mi ~ m^ 
such that there is no additional Am?, but only combinations of the new mixing angles. This 
case has recently been discussed by the MINOS collaboration [61]. Interestingly enough, in 
this limit we can still obtain a rather simple expression for P^^ with exact ^34 dependence, 
although with all other mixing angles set to zero and vanishing A21 [see App. A for details. 



^Notice that similar formulas, but using a different perturbative approach, have been derived in Ref. [53]. 



Eq. (23)]. From that formula, we can read off that there is, in principle, sensitivity to ^34 
via C34 in this limit, at least without correlations with the other mixing parameters. For the 
case IA41I ~ IA31I ^ A21, we have not found any useful analytical expression. We have to 
rely on numerical results. 

2.3 Simulation techniques 

We have performed the numerical simulation using the GLoBES software [63,64]. In order 
to use more than three neutrino flavors, we have defined a user-defined probability engine 
with the full four flavor probability including appropriate matter effects, assuming ?t,„ ~ ng. 
Instead of giving all details, let us just mention some of the complications. First of all, the 
GLoBES code (since version 3) is based on an efficient method to diagonalize Hermitian 
3x3 matrices developed by J. Kopp [65]. Since we had to diagonalize a 4x4 matrix now, we 
had to return to a GSL-based method, which is less efficient but easy to implement since the 
interfaces to GSL still exist from prior versions. Second, in order to avoid aliasing effects in 
the presence of fast oscillations, we have used the filter feature of GLoBES. And third, we 
have used the built-in standard minimization method GLB_MIN_NESTED_POWELL for 
the marginalization, which separates the systematics (using the pull method) and oscillation 
parameter marginalizations. In the presence of so many additional parameters, the new (but 
more efficient) method GLB_MIN_POWELL has rendered insufficient. An alternative may 
be the Markov chain interface MonteCUBES [66]. 

For our experimental setup, we mostly follow the International Design Study (IDS-NF) [18] 
baseline setup 1.0, which consists of 2.5 x 10^° useful muon decays per polarity, baseline 
and year with the parent muon energy Efj_ = 25 GeV. The total running time is assumed to 
be 10 years. Two magnetized iron calorimeters (fiducial mass 50 kton) are assumed at L = 
4000 km and L = 7500 km, respectively. The description of the Neutrino Factory is based 
on Refs. [18,25,67,68]. The considered oscillation channels are electron to muon neutrino 
(appearance channels) and muon to muon neutrino (disappearance channels) oscillations, 
unless noted otherwise. 

For the near detectors (we need two of them if the muons circulate in different directions 
in the storage ring [28]), we assume a fiducial mass of 32 1 and a distance oi d = 2 km from 
the end of the decay straight, which corresponds to the effective baseline of 2.28 km. This 
choice is motivated by a compromise between the sensitivity to as small Aml^ as possible, 
which prefers long baselines, and the longest possible baseline which seems to be plausible. 
At this distance, the peak sensitivity is achieved at around IAtti^j^I ~ 10 eV^ (c/.. Fig. 2 
in Ref. [54]). Note that the near detectors may, depending on the geometry and depth 
of the storage ring, also be built close to the surface if the up-going straight is used. We 
include electron and muon neutrino appearance and disappearance, which means that we, 
in principle, require charge identification. The muon neutrino channels are assumed to have 
at least as good capabilities as in the far detectors. The electron neutrino channels are 
assumed to have an efficiency of 40% and a charge mis-identification rate at the level of 1%, 
following Refs. [25,32]. Note, however, that the charge identification of electrons strongly 
depends on the detector technology and is, in general, difficult for high energies since the 
electrons produce quickly electromagnetic showers. From Eq. (11), however, we learn that 



Pee has a clean signature for 9^. This could be spoilt by V^e in Eq. (8) at the short 
baseline in the absence of charge identification if 624 was sufficiently large. However, this 
background is small for small mixings, and ^24 will be constrained by P^^, for which charge 
identification should be easier to achieve. Therefore, charge identification in the electron 
neutrino channels turns actually out to be not so important, as we have explicitly checked 
for the sensitivity limits by adding the electron neutrino and antineutrino event rates. As 
far as the beam geometry is concerned, at this distance the near detectors basically see the 
same spectrum as the far detectors and the effects of averaging over the decay straight are 
small at this distance [28]. 

In addition, we test the impact of optional OPERA-inspired (magnetized) Emulsion Cloud 
Chamber (ECC) i^T--detectors as near or far detectors in order to test the impact of the z/^ 
channels. The z/g — )■ Ur channel description is based on Refs. [25,69]. The z/^ — ?■ z/^ channel 
is assumed to have the same characteristics as in Ref. [69], governed by Refs. [36,53]. Since 
we assume that the hadronic decay channels of the r can be used as well, we assume a 
factor of five higher signal and background than in Ref. [69], i.e., 48% detection efficiency. 
In addition, we assume a detector mass of 10 kt, which is a quite aggressive choice to 
see where Ur detection may be really relevant. Note that is yet unclear if an ECC can 
be operated as close as 2 km to the Neutrino Factory because of the high scanning load. 
Therefore, alternative technologies may be preferable, such as a silicon vertex detector. In 
this case, other challenges have to be approached, such as the background from anti-neutrino 
charm production. These issues are currently under discussion within the IDS-NF. 

As far as systematics are concerned, any oscillation signals in the near detectors, especially 
in the electron and muon neutrino disappearance channels, are problematic in the presence 
of cross section uncertainties. Therefore, in Ref. [54], a set of near detectors at different 
baselines was proposed, similar to the Double CHOOZ or Daya Bay reactor experiments. 
With such a combination, cross section errors can be well controlled. Here we do not 
simulate such additional detectors explicitly, but only assume effective normalization errors 
of 2.5% (either from external measurements or from a near detector system). 

The best-fit oscillation parameters are taken as follows [70]: 

^12 = 34.4°, ^13 = 5.6°, ^23 = 42.8°, ^.4 = 0°; 

Amj^ = 7.59 X lO'^eV^ , Am^^ = 2.46 x lO'^eV^ , Am^^ = 1 eV^ (21) 

We impose external la errors on Am^i (4%) and ^12 (4%) and on Am|^ (10%) and ^23 
(10%) as conservative estimates for the current measurement errors [70]. We also include a 
5% matter density uncertainty [71,72]. 

3 Generalized exclusion limits 

In this section, we discuss general constraints to the new mixing angles 614, 624, and 634, 
and the additional mass squared difference Aml^ without any additional assumptions. In 
particular, we do not assume that Amli is in a particular range, such as the LSND-motivated 
one which leads to averaging at the long baselines. In addition, we do not assume that some 
of the not shown parameters take particular fixed values. 
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Figure 2: The exclusion limit for ain^ 29i4~Amli{i = 1,2,3) (region on r.h.s. of curves excluded). Here 
gives results of the standard IDS-NF (4000 km and 7500 km) setup with near detectors (thick solid curves), 
and the impact of the near detectors is also shown separately (thick dashed curves). The thin dashed curves 
only consider systematics. Mass ordering A assumed (c/.. Fig. 1), 90% CL (2 d.o.f). 



As performance indicator, we use the sensitivity to 6ij similar to the CHOOZ limit for the 
^i3-Ar?7,3j^ plane. We compute the simulated rates with 614 = 624 = 6^4 = and Am|j^ = 0, 
corresponding to the hypothesis of no effect of an additional sterile neutrino. Then we can 
discuss several exclusion limits for the new mixing angles, with the unknown parameters 
marginalized over. The exclusion limit for each new mixing angle will, in general, depend 
on Aml-^ similar to the CHOOZ limit. Therefore, we show three different exclusion planes 
di4-Aml^. Note that the sensitivity for any combination 9i4-9ki {i 7^ k) will typically vanish, 
since Aml-^^ is marginalized over. For a comparison to the existing literature, see the next 
section. 



Our main result can be found in Fig. 2: The exclusion limit for sin^ 29i4-Am 
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1,2,3) 



(region on r.h.s. of curves excluded). Let us first of all discuss the effect of the near detec- 
tors separately (thick dashed curves). Obviously, the main sensitivity is obtained at about 
Aml-^ ~ 10 eV^, which comes from the distance chosen for the near detectors. By changing 
the near detector locations, the position of the main peak can be controlled. If required 
by the recent MiniBooNE results, a longer baseline than 2 km may have to be chosen. As 
discussed in Sec. 2.2, the main sensitivity to 6*14 comes from electron neutrino (antineutrino) 
disappearance and the main sensitivity to 624 from muon neutrino (antineutrino) disappear- 
ance. Since the efficiencies for muon neutrino detection are typically better, the sensitivity 
to 624 is slightly better than that to 614 for our assumptions. As expected, there is no sensi- 
tivity to 634 coming from the near detectors, because the Ur disappearance channel does not 
exist. Note that the appearance channels at short baselines always depend on the product of 
two new parameters, c/.. Sec. 2.2. This means that no sensitivity can be obtained from the 
appearance channels without additional assumptions. Consider, for instance, V^^ in Eq. (8), 
which is sensitive to the product of ^14 and 6*24. In order to constrain sin^ 2^i4-Am4j^, ^24 
is to be marginalized over, which destroys the sensitivity. The problem is that this channel 
only limits a combination of 614, 624, and Am^^, which is difficult to display. Therefore, 
by our definition, we do not find any sensitivity from this channel, as opposed to Fig. 7 in 
Ref. [51]. 
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Table 1: Exclusion limits (90% CL, 1 d.o.f.) for several selected (fixed) values oi Arnli. Mass ordering A 
assumed. The simulation includes both near and far detectors. 

For the effect at the long baselines, it is first of all useful to consider the thin dashed curves 
with systematics only. In all three panels, the sensitivity changes as a function of Am^i 
in the region where Amli ~ Am"^^. It conies from the fact that the Neutrino Factory is 
sensitive to the atmospheric oscillation frequency, whereas for Aml^ ~ Am2i, no particular 
additional effects from the solar frequency can be found. As expected [cf., Eq. (18)], the 
main sensitivity is found for 6*24 which can be measured with the V^^ disappearance channel. 
However, there is also some sensitivity to 6*14, which vanishes after the marginalization, and 
some sensitivity to 6*34, which is even present for Am^^ = for systematics only [cf., the 
Appendix A for analytical formulas in that limit). After marginalization (thick solid curves), 
only the sensitivities to ^24 and ^34 remain in the Am^i regions close to the atmospheric 
Amli and above, where the effects of Am^i average out. Very interestingly, note that 
mixing angle correlations destroy the sensitivities for Amli = Amli, where m^ = m^ and no 
additional Am4^ is observable, leading to small gaps (see horizontal lines). The sensitivity 
to ^34 is not visible in Sec. 2.2 even for large Am\^. We have tested that it is a matter 
potential-driven, statistic limited higher order effect in ^34, present in the muon neutrino 
disappearance channels. Note that while one may expect some effect from V^r-, we have 
tested the impact of Vr detectors at all baselines, and we have not found any improvement 
of the sensitivities. The reason is that the ^34 effect at the long baseline comes with the 
same energy dependence as the 6*24 effect, which means that one cannot disentangle these, 
as discussed in Sec. 2.2. 

In view of the three panels, it is not easy to disentangle the parameters for arbitrary mas- 
sive sterile neutrinos. Parameter correlations lead to a pollution of the exclusion limit of a 
particular mixing angle with Am\]^. In addition, there is a competition between Am\]^ and 
Am|]^ at the long baseline. Near detectors, on the other hand, have very good sensitivities 
to 6*14 and 6*24 but cannot measure 6*34. Nevertheless, the absolute values of the sensitivi- 
ties are quite impressive, see Table 1. Especially, 6*24 can be very well constrained close to 
the atmospheric mass squared difference range. This indicates that sterile neutrino bounds 
in that range should be also obtainable from current atmospheric neutrino oscillation ex- 
periments. In addition, solar and supernova neutrino experiments may test even smaller 
Am\^ <C 10~^eV^, which is to be investigated. 

To illustrate the effect from different mass orderings, as shown in Fig. 1, we display these 
figures for the standard IDS-NF (4000 km and 7500 km) setup in combination with near 
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Figure 3: The exclusion limit for sin^ 26'i4-A77i4]^ (i = 1,2,3) (region on r.h.s. of curves excluded). Here 
gives results of the standard IDS-NF (4000 km and 7500 km) setup with detectors for the four different 
mass orderings in Fig. 1; 90% CL (2 d.o.f.). 

detectors in Fig. 3. Note tliat only tlie absolute value of the new mass squared difference 
is shown at the vertical axes. The upper peak hardly depends on the mass ordering, as it 
is obvious from the analytical short-baseline formulas in Sec. 2.2. The lower (long-baseline) 
peak, which is only present in the middle and right panels, somewhat depends on the mass 
ordering. We identify two qualitatively different cases: In schemes A and D, the sensitivity 
is destroyed just at the value of Am|]^. In these cases, c/.. Fig. 1, mass eigenstates 3 
and 4 are on top of each other, which means that there is no additional mass squared 
difference. The parameter correlations (marginalization over the unknown parameters) then 
destroy the sensitivity because the new mixing angles cannot be disentangled, in spite of the 
additional neutral current matter effect. This is different for schemes B and C, for which 
mass eigenstate 4 is on the opposite site of mass eigenstate 3. Although the absolute values 
of Amli and A?ti|^ are similar, these mass squared differences have different signs leading 
to different (charged current) matter effects. 

4 Exclusion limits with special assumptions 



Here we compare our analysis with different approaches in the literature under special 
assumptions. The most common assumption for sterile neutrino bounds at the long baselines 
is I A?724;^ I ~ 0{1) eV^, as motivated by LSND, which leads to averaging over the fast Amli at 
the long baselines. Another assumption, which we have found in Ref. [61], is Aml^ — t- 0. The 
third case we consider is the two- flavor short-baseline limit. These examples are particularly 
useful to discuss some subtleties when short- and long-baseline results are to be combined. 
Note that, for the comparison to the existing literature, we use studies which use the same 
parametrization of the four neutrino mixing matrix as ours for the sake of simplicity. 
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Figure 4: The exclusion limit for 6'24-6'34 in degrees (90% CL, 2 d.o.f.), where |Am|i| > |Ato3i|. We 
show the result for the IDS-NF setup (4000 km and 7500 km) using far detectors only (thick solid curve). 
The thick dashed curve corresponds to the result without correlations (systematics only), the dotted curve 
the result including a Vr detector at 4000 km. In addition, the current bounds are shown (from Ref. [60]). 



LSND-motivated Am^.^ 



The first assumption, wliicli we test, is an LSND-motivated iSm^^. In general, we assume 
^ |Am|]^|, whicli leads to tlie averaging of tfie fast oscillations at the long 



that 



|Am^i| 



baselines. This limit is frequently used for the discussion of long baselines only, z.e., without 
near detectors, since oscillation effects may be present in the near detectors. 

We show in Fig. 4 an example for such an analysis using the far detectors only (thick solid 
curve), which is to be compared with Fig. 9 in Ref. [53] for slightly different parameter val- 
ues (here the simulated ^14 = 0). For comparison, the current bound is shown (c/.. Fig. 2 in 
Ref. [60]). One can easily see that the Neutrino Factory could improve the current bounds 
on 6*24 and ^34 by a factor of a few if |A?Ti4]^| is assumed to be large. Note, however, that 
the marginalization over t^m^Y would lead to vanishing sensitivity. If the (thick dashed) 
curve with only systematics is considered, the cross terms in the long-baseline probabilities 
in Sec. 2.2, such as in "P^^ in Eq. (18), are switched off, and result improves along the 
diagonal. Such cross terms are also present in "P^^ in Eq. (19), sometimes called the "dis- 
covery channel" [53] , which leads to a slight improvement if an additional v^ detector at the 
intermediate baseline is used (thin dotted curve). However, as we noted above, the main 
effect on ^34 comes from the disappearance channels, especially of the very long (7500 km) 
baseline. 

Let us now test the impact of additional near detectors on this scenario, see Fig. 5. Ob- 
viously, we have to be careful, because a IS.m\^ ~ 1 eV^ might cause observable effects in 
the near detectors, and, in principle, we also have to marginalize over Aml^. We especially 
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Figure 5: The exclusion limit for ^24-^34 in degrees (90% CL, 2 d.o.f.)) where |Am|i| > |Ar 
Compared to Fig. 4, additional near detectors at 2 km are used. This leads to the dependence on the 
assumptions for Amli (and sin^ 26*14), as given in the plot legend (see main text for details). Here we assume 
sin^ 2^14 = 0.01 for the curves including near detectors. Note that the dashed-dotted curve corresponds to 
the thick solid curve in Fig. 4. 



expect some impact on tlie measurement of 6*24 by tlie muon neutrino disappearance channel 
in Eq. (14). If we fix Aml-^ = leV^, we find considerably better sensitivity including the 
near detectors, compare the dashed-dotted curve (without near detectors) with the thick 
dashed curves. However, note that the marginalization over Amli will in this case de- 
stroy the sensitivity again. This problem can be circumvented by additional assumptions. 
Consider, for instance, a non-zero value of sin^ 2^14 chosen by Nature. Then Aml-^ can be 
actually measured by Vee in the near detectors, see Eq. (11), and the marginalization over 
Amli can be performed. In this case, the impact of the Aml-^ marginalization is small, as it 
can be seen from the comparison between the thick dashed and thick solid curves. Another 
plausible assumption may be that Amli is so large that the effect even averages out in the 
near detectors. Therefore, we show the thin solid curve for a KeV sterile neutrino with 
Amli = 10^ eV^, which is sufficiently large. Here the near detectors somewhat improve 
the sensitivity compared to the no near detector case, but the effect is not as large as for 
the Aml-^ ~ leV^, sin^ 26*14 = 0.01 case. The purpose of this example is to illustrate that 
any combined fit of two new parameters not including Aml^, such as the one in Fig. 4, 
face subtleties if near and far detectors are combined. Typically, additional assumptions 
are needed, and the interpretation of the results becomes assumption-dependent. This is in 
contrast to the figures in the previous section, which do not depend on assumptions. 

A similar analysis is performed in Ref. [61] by the MINOS collaboration. The mixing matrix 
parametrization in this reference is equivalent to ours since they fix their 62, corresponding 
to our Si, to zero. One of the schemes tested in Ref. [61] assumed that |A?7i4i| ^ 
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Figure 6: Exclusion limits in the 034-^24 (a), 623-O24 (b), and ^34-^23 (c) planes at the 90% CL. The 
different curves correspond to MINOS and the Neutrino Factory (NF) with A?7i|]^ = 1 eV^. Here only one 
NF far detector at 4000 km is used (without near detectors). All contours represent 90% confidence level. 
The solid curves assume 613 = 0°, while the dashed curves assume ^13 = 12°. The best-fit values are marked 
in the figure, the one of 624 is zero. The MINOS curves are taken from Rcf [61]. 

(in mass ordering A). We sliow tlieir result for tlie combined fits in the 634-624, (a), 623-624 
(b), and 634-623 (c) planes in Fig. 6. Note that the not shown parameters are fixed, such as 
6*14 = and the phases, and that 6*13 is fixed to two different values (not marginalized over). 
The best-fit values are also marked. We also show the results for the Neutrino Factory 
under the same assumptions, where we use the 4000 km baseline only. Obviously, the 
Neutrino Factory would reduce the allowed parameter space significantly, especially if ^13 is 
large. Again, note that some of the parameters are fixed here, and the full marginalization 
would destroy the sensitivities. The impact of additional near detectors, which are especially 
sensitive to ^24, and an additional far detector at 7500 km, which is sensitive to 6*34, is shown 
in Fig. 7. Here the NF contours for large 613 = 12° almost shrink to points, and are hardly 
visible anymore. In this figure, Aml^ = 1 eV^ is assumed, and Aml^ is not marginalized 
over. For such a value of Amli, we would also expect a small effect in the MINOS near 
detector, which is, however, not considered in Ref. [61]. 



14 



NDFD@NF v.s MINOS 



15 



(a) 






60 r 

50- 

-n 

N 

40- 



-\-^ — I — I — I — \—¥ 



H — I — I — I — hM- 



(b) 




(c) 



Q3 



30 



40 50 

023 



MINOS, fl|V 
MINOS, fl|V 
NF, 613 = 12" 
NF, 611=0" 



20 40 

^34 



60 



Figure 7: Same as Fig. 6, but additional far detector at 7500 km and near detectors at 2 km included. 
The NF contours for large 6*13 = 12° almost shrink to points. 

The special case Am^.^ — >■ 

As similar kind of analysis can be performed in the special case Am|^ — )■ 0, for which also 
no additional mass squared difference appears and mi = m^. In this case, we have some 
sensitivity to ^34 in "P^/^, as shown in the Appendix A. The comparison between MINOS and 
the Neutrino Factory is shown for this special case in Fig. 8, and again the Neutrino Factory 
has an excellent sensitivity. The marginalization over Amli and the other parameters will, 
however, destroy the sensitivity to 6*34, since there is no sensitivity in Fig. 2 for Am|^ — >■ 0. 
This is another good example that the set of assumptions determines the outcome, whereas 
the general analysis would simply produce no sensitivity. 

Comparison to two-flavor limits 

Especially for short-baseline new physics searches, in the literature often two baseline limits 
are discussed. That is, at a short baseline where no standard oscillations occur, sterile 
neutrinos are searched for using a two flavor fit. An example is the short-baseline electron 
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Figure 8: Exclusion limits in the 934-623 plane at the 90% CL. The different curves correspond to MINOS 
and the Neutrino Factory (NF) with Amli = eV^. Here only one NF far detector at 4000 km is used 
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and 7500 km together with the near detectors at 2 km are included. The solid curves assume 613 = 0°, 
while the dashed curves assume 613 = 12°. The best-fit values are marked in the figure, the one of 6*24 is 
zero. The MINOS curves are taken from Ref. [61]. 



neutrino disappearance considered in the two flavor limit: 

., , ^, ,, /Am'^L 
sin^(2e) sin^ 



AE 



(22) 



In Ref. [54] , even the test of CPT invariance at such a short baseline was considered. From 
our discussion in Sec. 2.2, it is clear that the short-baseline disappearance channels are the 
key measurements for the active-sterile-mixings, because they are uniquely sensitive to one 
parameter (or one mixing matrix element). On the other hand, the appearance channels, 
such as the tau neutrino appearance in NOMAD and CHORUS or the electron neutrino 
appearance in MiniBooNE, measure combinations of parameters. A signal in this channel 
would inevitably point towards new physics, but these channels are not very useful to limit 
individual active-sterile mixing parameters (or mixing matrix elements). 

Eq. (22) can be directly related to our parametrization with the identiflcation 6 — )■ 614 
and Am^ — )■ Amli. Correspondingly, the results of our analysis for ^14 are similar to 
Ref. [54] (see Fig. 3). On the other hand, it was emphasized in Ref. [54] that cross section 
measurements will spoil the disappearance measurements, and that these can be controlled 
with additional near detectors very close to the source. Since there is no CP violation in 
short-baseline disappearance measurements (neither from CP phases nor from the MSW 
effect), the disappearance channels can also be used for clean CPT invariance tests. 

Note that a CPT invariance test of the MiniBooNE result would be, in principle, also 
possible at the Neutrino Factory, because both P^g and P^e could be tested. However, 
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for these appearance channels electron charge identification is required to distinguish, for 
example, Pp^g from the Pee background. Therefore, the direct test of the recent anomaly [11] 
may be difficult. 

5 Summary and conclusions 

In this work, we have discussed sterile neutrinos beyond LSND, i.e., sterile neutrinos with 
small active-sterile mixings and an arbitrary Amli. We have used the Neutrino Factory for 
the simulation, since this experiment can be used for a self-consistent approach including 
near and far detectors. We have used the simplest possible hypothesis, namely one extra 
sterile neutrino in a 3+1-like scheme in which the standard scenario is recovered for small 
active-sterile mixing angles. While we obtain the expected sensitivity to the active-sterile 
mixing in the large LSND-motivated |A?72|^| > 1 eV^ region, we also find sensitivity close 
to the atmospheric Amg^. We have pointed out that there are no global fits for this case 
yet, and a re- analysis of atmospheric and solar data should demonstrate what we can learn 
for Amli <^ 1 eV^. Note that recent cosmo logical fits point towards one or two such light 
sterile neutrinos [12]. 

We have demonstrated that, especially at the short baselines, the disappearance channels 
are the primary channels of interest for light sterile neutrino constraints, no matter if a 
parametrization-independent approach or a particular parametrization for the mixing ma- 
trix is used. For the combined analysis of short and long baselines including charged and 
neutral current matter effects, however, we have used a particular parametrization. We have 
demonstrated that the most general constraint on sterile species can be shown as exclusion 
limits in the ^j4-A?n4^-planes (z = 1, 2, 3), similar to the CHOOZ limit for Oi^-Amli. The 
Neutrino Factory turns out to have excellent sensitivity to the three mixing angles in a wide 
range of Aml-^. However, one of the three mixing angles (^14 in our parametrization) can 
only be very well measured for large |A?Ti4^| > 1 eV^ at the near detectors, and one {634^ 
in our parametrization) better for small |A?7i4]^| ~ |Am|]^|, whereas the third (^24) can be 
measured in the combined range. An electron neutrino disappearance channel at the long 
baselines could solve this problem for 614, which is, however, difficult at the Neutrino Fac- 
tory, because electron charge identification might be required. An improved measurement 
of 634 for large |A?Ti4^| > leV^ would require a hypothetical Ut disappearance channel. We 
have also investigated the impact of different mass orderings (c/.. Fig. 1) on the sensitivities. 
We have found that there is a qualitative difference between the cases Am^i ~ Am|^ and 
Am|^ ~ —Amli due to different matter effects. Furthermore, we have tested the impact of 
additional Uj. detectors with an aggressive 10 kt OPERA-like detector at the short and long 
baselines, and we have not found any significant effect on the sensitivities. 

Apart from the general constraints, we have compared our analysis to special cases in the 
literature. For instance, we have tested the case | Am4]^| > 1 eV^ leading to averaging of the 
fast oscillations at the long baselines. We have looked into the special case of a 62A-O34: fit. We 
have shown that this assumption involves subtleties, especially if additional near detectors 
are considered, because the assumption for Aml-^ implies that there may be effects in the 
near detectors. We have illustrated that in this case additional assumptions are required. 
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and the interpretation of the sensitivities becomes strongly assumption-dependent. This 
does not apply to our general bounds we discussed above. We have also tested the impact 
of Vt- detection ("discovery channel"), and we could find a marginal improvement. The main 
sensitivity, however, comes from the muon neutrino disappearance channel for E^ = 25 GeV. 
We have moreover compared the Neutrino Factory to a recent MINOS analysis using the 
same assumptions, and we have found excellent sensitivities. However, including the full 
marginalization over the unknown parameters, the sensitivities were basically destroyed. 

In conclusion, although the sensitivities of the Neutrino Factory to light sterile neutrinos 
in the atmospheric mass squared range may not be extremely compelling, these bounds 
are very robust and independent of any special assumptions. If special assumptions are 
used, such as often done in the literature, the sensitivities look extremely good. Of course 
our results are not the most general case, one could also study 3 + A^ scenarios. However, 
even for the 3 + 1 scenario, this analysis is extremely challenging because there can be 
interference between the known and the unknown mass squared differences. A re-analysis 
of atmospheric and solar experiments could provide clues on what one can learn already 
from existing experiments. Finally, even the test in supernovae may be interesting, because 
one may expect additional resonances or new (neutral current) matter effects. 
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A P^^ in the limit IA41I < A31 

In this appendix we show the transition probability P^^ in the limit IA41I ^ A31 and 
A21 = 0. We study the explicit dependence on 6*34 to all orders in perturbation theory, 
while keeping 623 = 45° and all other mixing angles vanishing. We get: 



'^'^ C2 (1 + cl^y 



2C34 sin C(A„ sin A„ — C34A31 sin A, 
7TT\ ', J2T9. 



(23) 



C(l + ci^2 



^34 J 



18 



where 

C = ^24,A,,A^ + Al, + A2 . (24) 

We see that the dependence on 634 (through the C34 terms) appears in terms with different 
energy dependencies, which malces a spectral analysis at the Neutrino Factory useful to 
get a non- vanishing sensitivity even in the IA41I -C A31 limit. Interestingly enough, the 
matter effects are driven by the neutral current potential A„, because of the vanishing 6*13 
approximation . 
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